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Abstract: Open questions remain after the confirmation of LMA as the dominant solution 
to the solar neutrino deficit. These are the apparent time modulation of the solar neutrino 
event rates in the Homestake, Gallium and SuperKamiokande experiments, possibly re- 
lated to solar magnetic activity, the discrepancy between the event rate in the Homestake 
experiment and its theoretical prediction and the absence of the electron spectrum upturn 
in SuperKamiokande at energies below 6-8 MeV. We search for a possible understanding o 
of these questions in the context of resonant spin flavour precession to sterile neutrinos, 
assuming a class of magnetic field profiles anchored in the upper radiation/lower convection 
zone. We consider the simplest such model beyond the standard 2v flavour LMA one, with 
one single magnetic moment transition between active and sterile state and vanishing vac- 
uum mixing. The preferred mass square difference is Amf = 0(W~ 8 eV 2 ). The Ga rate 
appears to be the most sensitive of all to solar activity. It is also found that a field profile 
extending within a longer region in the radial direction is favoured over another with a 
shorter span, and leads to a stronger suppression than in the LMA case of the intermediate 
energy neutrinos and some of the S B ones. 
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1. Introduction 

The KamLAND experiment [|| seems to confirm the Large Mixing Angle (LMA) solu- 
tion to the solar neutrino deficit, originally presented in 1992 [g], and to rule out spin 
flavour precession (SFP) ||. Even the possibility of a sub-dominant SFP is substantially 
constrained j| in view of the recent upper bound on solar antineutrinos corresponding to 
0.028% of the S B neutrino flux obtained by the KamLAND collaboration ||. It would be 
however premature to assert that LMA || is the final and complete solution, at a time 
when solar neutrino experiments enter a new era of precision measurements and there are 
hints of possible phenomena which fail to fit in this scenario. 

One of the generic predictions of LMA is an exceedingly large rate |7j for the CI 
experiment [pi]. This discrepancy reaches an excess of 2.5 a if one uses the latest BP'04 



results |J. Moreover, despite accurate predictions for the SuperKamiokande [10] and SNO 



11|] total reduced rates, the corresponding LMA spectrum predictions exhibit an upturn 
below (6-8) MeV which is absent in both data. While the absolute values of the rates 
may not be significant, because they involve a normalization to a specific 8 B flux whose 
value is known within a (15-20)% theoretical uncertainty, the absence of this upturn in the 
data may motivate new physics beyond LMA, since it does not rely on normalization to a 
specific model 1 . Reconciling the theory with the data in this aspect implies a decrease in 
the survival probability for neutrinos with energies in the range (0.8-8) MeV, that is, those 
which are observed by Chlorine and unobserved by SK and SNO. 

Another hint which may point to new physics lies in the possible time dependence 
of the neutrino signal. Solar neutrino statistical analyses confronting theoretical model 
predictions with data have been using time averaged event rates or fluxes || |l2| and could, 
for this reason, be missing important information. The examination of the data on a time 
basis by the Stanford Group has in fact provided increasing evidence that the solar neutrino 



flux is not constant, but varies with well-known solar rotation periods |13|, 14, 15]. Such a 



1 The present statistics both in SK and SNO are however not enough to perform a strict statement 
concerning acceptance or rejection of a model only on grounds of this upturn. 



situation, if confirmed, can in no way be explained by the LMA scenario. It can neither 
be understood in terms of an SFP transition of v e into or v T Q , as this would originate 
a sizable and unseen v e flux. 

The intermediate energy dip in the survival probability and the possibility of modula- 
tion of the neutrino signal constitute our main motivation to investigate scenarios with a 
sterile neutrino. A magnetic moment driven conversion from active to sterile neutrinos and 
the fact that these leave no trace in the detectors appears to be an attractive possibility to 
generate these two effects. In fact, time dependence requires this conversion to be related 
to solar magnetic activity or to solar rotation with an axially asymmetric magnetic field 
profile. 

Our purpose in this paper is to perform a combined prediction of all solar neutrino data, 
namely the rates for Chlorine |8| and Gallium []lq], [17| experiments, the SuperKamiokande 
p0| ] reduced rate and spectrum and the reduced rates in SNO data p| . We will not be 
concerned at this preliminary stage with x 2 fittings, but rather with showing the change in 
the LMA probability shape and solar neutrino predictions resulting from a possible time 
dependent magnetic field related conversion to the sterile neutrino. In section 2 we present 
the general structure of the model, in section 3 we derive the numerical predictions of all 
rates and discuss other consequences of the model for two specific solar field profiles and 
in section 4 we draw our main conclusions. 



2. Structure of the Model 



In order to expound our model, we consider a system of two active neutrinos and a sterile 
one which mix in the mass eigenstates i>o, v\ and V2- This can be parametrized by the 
following rotation matrix [18| 
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with 6 denoting the usual LMA vacuum mixing angle and a the vacuum sterile mixing. 
A straightforward but tedious calculation then leads to the following form of the vacuum 
Hamiltonian 
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with obvious notation. We seek for a model which should be the simplest departure from 
the conventional LMA case and be able to generate a time dependent transition into a 
sterile neutrino. To this end we take a vanishing mixing angle, a=0, so that active states 
v e , communicate to the sterile one through magnetic moment transitions only. The 



matter Hamiltonian is thus 
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with V e and being the matter induced potentials for v e and respectively and fii, ^2 
their transition moments to the sterile neutrino. We will consider field profiles concentrated 
around the bottom of the convective zone, as motivated by the dynamo theory 14, 15], 



and hence require spin flavour precession to be resonant in this region. Therefore the two 
processes (LMA and RSFP) occur sequentially at very different solar radii, with the RSFP 
critical density determined by a mass square difference between one of the active states 
and the sterile one 0(W~ 8 eV 2 ). An efficient conversion to the sterile neutrino requires 
that the transition moment associated with this order of magnitude difference dominates 
over the other. Given these facts we choose Amf = 0(W~ 8 eV 2 ) and so ^2 = 0. 2 

We obtain the neutrino survival probability {P(y e v e )=P ee ) by assuming the prop- 
agation inside the sun to be adiabatic except in the region where the adiabaticity pa- 
rameter reaches its minimum. This is essentially a generalization of the (2x2) case pro- 



posed some time ago [19|. In fig.l we show the evolution of the mass matter eigenval- 
ues inside the sun for a typical case. It corresponds to a neutrino energy E = 2.1MeV, 
Amf = -1.7x 10 -8 eV 2 and field profile 2 described below with peak field B = 3.0 x 10 5 G. 
Of the three eigenvalues, it is seen that v\ and vq become the closest in the convection zone 
{r/Rs — 0.8), so that adiabaticity may be broken in the transition at this critical density. 
We recall that a strong field in this region will provide an efficient, adiabatic conversion 
corresponding to a larger eigenvalue separation while a weaker field will generate less con- 
version, with closer eigenvalues. Moreover the LMA oscillation and the v 2 — ¥ ^0 transition 
are strongly adiabatic. This suggests that a field profile exhibiting a time dependent peak 
value deep inside the convective zone will manifest itself in the modulation of the neutrino 
signal. Hence we need only to consider the jump probability Pq between v\ and and 
the expression describing the three neutrino propagation is 
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Figure 1: The evolution of the mass matter eigenvalues along the solar neutrino tra- 
jectory. The example chosen is for profile 2 (see fig. 2) with neutrino energy E = 2.1MeV, 
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no 



= -1.7 x 10~ H eV 2 



best fit values JTl[: Am|i 
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Here Hp denotes the diagonalized Hamiltonian, r$ the neutrino production point, r D the 
detection point and r x the location of the RSFP critical density. [R] is the rotation matrix 
from the weak eigenstate to the matter eigenstate basis. Its form in the vacuum is given 
by the transpose of (2.1) with a = 0. In matter with magnetic field it becomes in its most 
general form 
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where 6q\, 602 ar e the mixing angles induced by the neutrino magnetic moments fii, ^2 
respectively. The matter mixing angle 9\2 becomes in the vacuum the LMA mixing angle 8, 
eq. (2.1). In equation (2.4) [Ri\ converts the weak eigenstate v e produced in the solar core 
into a combination of Hamiltonian eigenstates which propagate adiabatically (matrices [C] 
and [A]) except for the correction at the critical density described by the matrix [B] where 



2 The other situation, namely a non- vanishing fi2 and /ii = would correspond very closely to the LMA 
one and two sizable transition moments to an intermediate possibility. 
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u\ and vq nearly meet. In order to account for level crossing at the resonance between 
states z^o, vi, the corresponding first two rows of the rotation matrix R are interchanged 
so as to obtain R Q . Hamiltonian eigenstates are then converted back to the weak basis 



through Rj at the detection point. Taking #02 
[Ri] and [R%] therefore read 



(/j, 2 = 0) and using eq.(2.8), matrices 
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For the jump probability Pq we use the Landau Zener approximation 
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(2.11) 



with the adiabaticity parameter given by the ratio between the corresponding eigenvalue 
difference and the spatial rate of the mixing angle 
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Using equations (2.4) - (2.10) and the fact that at the neutrino production and detection 
points the absence of magnetic field implies the vanishing of angle 9qi, the expressions for 
the survival and conversion probabilities become 
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(2.13) 
(2.14) 



where we neglected fast oscillating terms which average to zero and subscript % refers to 
the mixing angle at the production point. The LMA survival and conversion probabilities 
in the 2-flavour case can be readily obtained from eqs. (2.13), (2.14) for Pq = 1 (vanishing 
magnetic field at the resonance point). 

The expression for the SK spectrum to be used is 
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with standard notation J20| and a similar expression for the SNO one. We take the SK and 



SNO energy resolution functions from |21|,|p2| 
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Figure 2: Solar magnetic field profiles. These are motivated by the dynamo process near the 
tachocline J]l| : they exhibit a common and fast rise starting at the upper radiation zone reaching 
a peak at the bottom of the convection zone. Profile 2 extends over a longer region, hence it leads 
to more efficient sterile neutrino conversion. 



3. Solutions and Discussion 

The two classes of field profiles to be used are characterized by a sharp increase starting in 
the radiative zone just below the tachocline, reaching a peak at the bottom of the convective 
zone and a smoother decrease up to the solar surface. They are (see fig. 2) 
Profile 1 

B = ZZZUT^TZ ITTT ' x c < x < 1 ( 3 - 2 ) 



cosh(18(x — x c )) 



Profile 2 



B = B °\ 1 -{^jt) I ' x n< x < x c (3-3) 



B = B Q \l-[ Y — 2- ) ) , x c < x < 1 (3.4) 



with vanishing field for x < x R . Here x is the fraction of the solar radius, x R = 0.65, 
x c = 0.71 and we take in all cases a magnetic moment fi\ = 10~ 12 /ib- The preferred peak 
field values are found to be in the range Bq = (2 — 3) x 10 5 G. 
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Figure 3: Survival probabilities and SuperKamiokande spectrum for profile 1 (curves 

1 and 2). For comparison, the LMA curves are also given (curves 3). The dip in the probability 
hardly affects the shape of the spectrum. (See also table II). 



In table I we show the solar neutrino data on rates and their comparison with the SSM 
predictions including the theoretical uncertainties 3 . We refer to the parameter values 
as given from global fits [11| Am^ 



7.1 ±5:§xio- 



^eV 2 , 9 = 32.5±2'| degrees. Our main 



results are summarized in figs. 3 - 6 and tables II, III. 



Experiment 



Data 



Theory Data/Theory Reference 



Homestake 
SAGE 

Gallex+GNO 
SuperKamiokande 
SNO CC 
SNO ES 
SNO NC 



2.56 ±0.16 ±0.15 
70.9 ±|;1 ±|; 2 7 
70.8 ±4.5 ±3.8 
2.35 ±0.02 ±0.08 

i cq _t_0.08 _t_0.06 
l.oy 3=0.07 ^O.OS 

2.21 ±°;jg ±0.10 
5.21 ±0.27 ±0.38 
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5.82 ± 1.34 



0.301 ±0.069 

n ^41 +0.070 
u.o4in= 062 

n c; 4n+ 0.067 
u - O4U=t 0.061 

0.404 ± 0.094 
0.273 ± 0.065 

0.380±8i83 
0.895 ±0.221 
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Table I - Data from the solar neutrino experiments. Units are SNU for Home- 
stake and Gallium and 10 6 cm _2 s _1 for SuperKamiokande and SNO. The uncertainties in 
the fourth column include solar standard model errors. 



3 Here we use the BP04 standard solar model 
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Figure 4: Same as fig. 3 for profile 2. A downturn of the spectrum appears in this case as some 
S B neutrinos (E > 5MeV) are converted to sterile ones. (See also table III). 
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Table II - Profile 1: Model predictions for rates. Units are SNU for CI, Ga and 

the ratio of the model prediction by the SSM value for SK and SNO. See also fig. 3 for the 
survival probabilities and SK spectrum. Here Am^ (eU 2 ) = 6.5 x 10~ 5 eU 2 , 9 = 30.9°. 

Fig. 3 (profile 1) and fig. 4 (profile 2) show the electron neutrino survival probability 
and the electron energy spectrum in SK for two values of Am 2 in each profile together 
with the LMA probability and spectrum, for comparison, near the pure LMA best fit point. 
It is seen that the upturn in the LMA spectrum can be reduced or made to disappear. A 



flat spectrum, as suggested by the data [ 10 1 , is in fact obtained for a convenient choice 



of Am 2 and Bq in the case of profile 2 but not for profile 1. For both profiles, however, 
the CI rate is substantially decreased from its LMA prediction of approximately 3.1 SNU 
(see tables II, III). Parameter Am 2 also determines the location of the important ui — > uq 
resonance, with less negative values corresponding to resonances closer to the solar surface. 
Also the longer spatial extension of a strong field in the case of profile 2 implies more 
conversion efficiency, a fact which is clearly reflected in a wider dip in the probability curve 
and even a flatness in the spectrum for Am 2 =— (1.6 — 1.7) x 10 -8 eV 2 . All rates can be 
brought to within 1<t of the experimental value with theoretical errors included. In view 
of our motivation to account for possible modulations in the rates, we will not attempt at 
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quantitative and detailed fits at this stage, since the solar neutrino data used are averages 
and further, they do not refer to the same periods in different experiments. 
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Table III - Same as table II for profile 2. See also fig. 4 for the survival probabilities 
and SK spectrum. Here Am^ (eV 2 ) = 6.8 x 10" 5 eV 2 , 9 = 32.1°. 

We note that the predictions from profiles 1 and 2 exhibit a manifest difference: in 
fact for profile 2 increasing the dip in the probability through an increase in the peak field 
leads to a stronger downturn of the spectrum, while for profile 1 an increase in the peak 
field leaves the spectrum unaffected (figs. 3 and 4). This is because in profile 2 a higher 
peak field leads to a stronger suppression of the low energy sector of the 8 B neutrinos, 
which are present in the SK rate, along with pep and CNO ones. For profile 1, which is 
more localized, an increase in the peak field affects only the suppression of the intermediate 
energy neutrinos in CI, especially the pep neutrinos with energy E pep = 1.4AM eV and the 
CNO ones, all absent in the SK rate and accounting for only a minor fraction of the CI 
rate. 

In figs. 5, 6 we show the dependence of rates on the peak field value for profiles 1, 2, 
from a starting point with unaffected LMA predictions, up to Bq = 3 x 10 5 G (2^] with a 
view on a possible anticipation of future studies of time dependence in neutrino signals. 
Four different cases of Am 2 in the relevant range (10~ 8 eV 2 ) are considered (see table IV). 
The values of Am§i and 9 are slightly different in figs. 5 and 6, hence the small discrepancy 



in the rates in the pure LMA limit. They are however well within their la range [11] 
(see the captions of tables II, III). We find that the Ga rate is the most sensitive of all to 
changes in Bq with a maximum modulation for profiles 1 and 2 respectively 

^Ga „ 22% ARca „ 4Q% (3 5) 

Rgcl Rgcl 

At this point it should be noted that in their detailed analysis, the Stanford Group [jl3| finds 
the Ga rate to exhibit the most manifest time dependence of all with two clear peaks in the 
ranges 45-75 and 90-120 SNU. It is unclear at this stage whether this apparent convergence 
of results is essential or not, a situation which obviously deserves further investigation. 





a 




b 




c 




d 


Profile 1 


-0.5 x 10" 


- 8 eV 2 


-1.35 x 10" 


S e y2 


-1.65 x 10" 


S e y2 


-5 x 10~ 8 eV 2 


Profile 2 


-0.5 x 10" 


-s eV 2 


-1.6 x 10" 


'eV 2 


-1.7 x 10" 


8 e y2 


-5 x 10~ 8 eV 2 



Table IV - The values of Am 2 labeled a, b, c, d in figures 5 and 6. 
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Figure 5: Event rates as a function of the peak field Bq for profile 1. Variations of solar 
activity, manifested in a variable field intensity, are reflected in variations of the event rates. SK, 
SNO-CC and SNO-NC panels are all within la of the expeimental value including theoretical errors. 
The LMA parameters are slightly shifted from their best fit values 111). (See also table IV). 
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The CI rate is less sensitive with a maximum 

^ ~ 12% , ^ ~ 17% (3.6) 
Hci tici 

for profile 1 and profile 2 respectively. For all other rates this sensitivity is always below 
7% for profile 1 and 15% for profile 2. All these modulations depend of course on the value 
of Amf chosen. Furthermore the relative sensitivities of different rates to Am^ change 
from one rate to another. The cases described in detail in tables II, III and figs. 3, 4 can 
also be seen in a different context in figs. 5, 6. 

In all our calculations we have used a neutrino magnetic moment \i v = 1CP 12 ^£- This 



is larger than the supernova bound [24, 25] by an approximate factor of 2 which applies 



for active-sterile magnetic moment conversions. Hence, if strictly enforced, this bound 
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Figure 6: Same as fig. 5 for profile 2. The LMA parameters are slightly shifted from their best 
fit values (See also table IV). 



would require a peak field value of (4 - 6) x 10 5 G instead of (2 - 3) x 10 5 G. However the 
supernova bound is plagued by uncertainties in the models of the supernova core p5| or 



based on rough estimates of the supernova energetic [24], so it cannot be expected to be 
quite stringent. On the other hand the solar magnetic field profile is poorly known, both 
in shape and strength. Whereas a field of the order of 3 x 10 5 G is quite possible in the 
tachocline ] [23|] , it is not clear whether this can be exceeded by an extra factor of (1.5-2) or 
whether this factor could come from intrinsic inaccuracy of the supernova bound. 



4. Conclusions 

To conclude, now that KamLAND has established that LMA is the dominant solution 
to the solar neutrino problem, there remain a few questions, which may turn out to be 
quite important, as LMA may be incomplete and new physics might be necessary. To this 
end, the hints from solar neutrinos are a possible time variation of the neutrino event rate 
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1 13, 14, 15] in the experiments, the absence of an electron spectrum upturn for low energies 
in opposition to the LMA expectation, and a rate in the CI experiment lower than the 
LMA predicted one. 

We proposed an answer to these questions by adding a sterile state to the two solar 
neutrino system. Such scenarios have already been developed in the literature [18, 26 1 in 
the context of oscillations alone, in particular it was shown that a sterile state may provide 
a solution to the above questions [18|, except for the possible time modulation. Visible 
states communicate in those models with the sterile one via a vacuum mixing angle. In the 
model proposed in this paper, they communicate instead with the sterile one via magnetic 
moment transitions, so that time dependence of neutrino signals, if confirmed, can be 
directly connected to solar magnetic activity or to solar rotation. The new mass square 
difference is about three orders of magnitude below the LMA one, so that the transition 
to the sterile state resonates mainly in the tachocline, where the solar magnetic field is 
expected to be concentrated. In this way, a situation in which the transition moment 
associated with the smallest mass square difference dominates over the other leads to the 
most efficient conversion. 

The upturn in the spectrum can remain unchanged, be reduced or eliminated, depend- 
ing on the choice of this mass square difference and the field profile. As far as the CI 
rate prediction is concerned, a sizable reduction from its LMA value is always obtained 
in the class of field profiles investigated. Profiles 1 and 2 illustrate these situations: the 
elimination of the spectrum upturn may or may not accompany the reduction in the CI 
rate (profiles 2 and 1 respectively). It should be noted however that both SK and SNO 
data do not contain at present enough statistics that allow a strict statement on the exact 
spectrum shape. 

We have also investigated the evolution of the event rates with the peak field values in 
each of the two profiles and found the Ga one to be the most sensitive of all. Interestingly 
enough, it was found by the Stanford Group from a combined analysis of Gallex-GNO 
and SAGE that the Ga rate exhibits a strong time dependence with two clear peaks at 
45-75 and 90-120 SNU. This apparent convergence between ours and their results certainly 
deserves further investigation. The CI rate comes next to the Ga one in sensitivity, followed 
by all others (see figs. 5, 6). In all cases we used the convenient mass square difference order 
of magnitude, Amf = O(l0~ 8 eV 2 ). Profile 2, with a longer spatial extension, leads to 
more sensitivity of the rates. We have not performed x 2 fittings, as, within the motivation 
of the present paper, they will be justified only for time dependent solar neutrino data 
which may soon become available. 

Finally, the uncertainties both in the supernova neutrino magnetic moment bound 
and in the solar field profile make it premature to discard or rule out the present analysis. 
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